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ABSTRACT 



Spectropolarimetry of Broad Absorption Line Quasars (BALQSO's) has 
demonstrated that the geometry of the absorbing material is far from spherically 
symmetric. Calculations of accretion disk spectra suggest that the intrinsic radiation 
pattern of quasars is also anisotropic. Because the quasar counts distribution is very 
steep, if the orientations of these two anisotropies are correlated, optical flux-limited 
samples would be very strongly biassed with respect to discovery of BALQSO's. In 
particular, currently favored models suggest that BALQSO's may be much more 
common than in the samples. If so, the intrinsic covering fraction of absorbing matter 
must be relatively large. Numerous consequences follow, including prediction of a new 
population of hard X-ray sources (as may have already been seen by ASCA), and 
a possible explanation of the anomalously strong NV 1240 lines seen in many high 
redshift quasars. 



Subject headings: accretion, accretion disks — galaxies:active — quasars: absorption 
lines — quasars: emission lines — quasars:general 
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1. Asymmetric Absorption and Anisotropic Intrinsic Radiation 

Only about 10% of quasars found in optical flux-limited samples exhibit broad absorption 
lines, resonance line absorption troughs extending ~ 0.1c to the blue of the emission line 
centers (Weymann et al. 1991). Because in most other respects Broad Absorption Line Quasars 
(BALQSO's) are quite similar to ordinary quasars, it is generally thought that rapidly-moving 
absorbing matter exists in all quasars, but occupies only a fraction of the solid angle around them 
(Weymann et al. 1991). The fact that the remaining light in the troughs is often strongly polarized 
indicates that the obscuration is not randomly distributed around the quasar in a statistically 
spherically symmetric manner, but is instead strongly aspherical. Equatorial wedges are the most 
popular geometry used to explain the observations (Goodrich & Miller 1995; Cohen et al. 1995; 
Brotherton et al. 1997). 

Quasars are also likely to be anisotropic in their intrinsic continuum emission. The most 
plausible source for the photons they produce in the rest-frame optical and ultraviolet bands is 
an accretion disk, although detailed models for this emission are not altogether successful (see 
Blaes 1998 and Krolik 1998 for reviews). While the specific character of the intrinsic angular 
radiation pattern is highly model-dependent, in virtually every model the anisotropy is significant. 
Typically, the wavelength range around 1000 - 2000A is limb-darkened, while shorter wavelengths, 
produced closer to the black hole, are limb-brightened (Sun & Malkan 1989, Laor & Netzer 1989, 
Agol 1997). 

If the orientations of the intrinsic radiation pattern and the absorbing matter are correlated, 
optical flux-limited samples will be biassed with respect to the detection of BALQSO's (cf. 
Goodrich 1997, who suggested that there is bias due to extra attenuation on the lines of sight 
through the absorbing matter). This effect is strengthened by the steepness of the quasar counts 
distribution. 

For example, suppose that the differential source count distribution if quasars radiated 
isotropically were dN/dF* = KF~ a , where N(> F*) is the number of quasars per sky solid angle 
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with flux greater than isotropic flux . Then the distribution per cosine of the viewing angle is 
d 2 N/dF*dn = (K/2)F~ a . If we now allow for anisotropy by setting F j, s = F*f(n) (implicitly 
assuming that the angular radiation pattern is independent of all other quasar characteristics such 
as luminosity or redshift), the differential source count distribution that we would actually observe 
is 

9N f 1 , 8 2 N yi d 2 N dF* K p x 

For optical fluxes in the neighborhood of = 16 - 17, a ~ 3.25 (Hartwick & Schade 1990). 
Since a quasar is only called a BALQSO if it is found within an equatorial wedge defined by 
\fjb\ < //*, the BAL fraction in a flux-limited sample is 

= jy dfid 2 N/dF obs dfi = !t djxT^i} 
IBAL Jo 1 dfid 2 N/dF obs dfi Jo 1 d^f-Hp) ■ 

Thus, if quasars are dimmer in the BAL direction, their population fraction at fixed observed flux 

is diminished because they are submerged in the larger population of intrinsically fainter non-BAL 

quasars. 

As already remarked, we cannot compute /(//) with great confidence. However, two plausible 
guesses will illustrate the possible magnitude of its effect. The first guess is the simplest possible 
model: that the local limb-darkening is given by the form appropriate to a gray atmosphere in the 
Eddington approximation (Mihalas 1978), but that the integrated flux is further multiplied by a 
factor of jJL because the radiation comes from a flat disk. Then = + 1.5|//|). As a foil to 

this form, we use a rough analytic fit to detailed radiation transfer calculations which also include 
a complete treatment of general relativistic corrections (E. Agol, private communication). This 
analytic fit is / 2 (M) - 0.23 + 1.54|/x|. 

The predicted BAL fractions are shown in Fig. 1. Generically, any significant limb-darkening 
in the direction of the absorbing matter leads to a very strong bias against the discovery of 
BALQSO's. In order to find ~ 10% BALQSO's in an optical flux-limited sample of limb-darkened 
sources, the true covering fraction /i* ~ 0.5. When /U* exceeds ~ 0.5, the fraction of BALQSO's 
discovered in a flux-limited sample rapidly approaches unity as the limb-darkening weakens and 
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the fraction of non-BALQSO's falls. Conversely, if the absorbing matter is correlated with the 
bright direction, the true fraction might be substantially less than the fraction inferred from 
flux-limited samples. 




Fig. 1. — The fraction of BALQSO's that would apear in an 
optical flux-limited sample for three different limb-darkening 
laws. The solid curve is the prediction of gray atmosphere 
limb-darkening on a disk, the dotted curve the approximation 
to detailed disk models, and the dashed line gives Jbal if the 
continuum were radiated isotropically. 



This picture also has one more corollary regarding the BAL fraction: if the absorbing matter 
is in the dim direction, the very faintest quasars (i.e. those at flux levels where a falls below unity) 
should be disproportionately BALQSO's. 
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2. Consequences of Large Covering Fractions 

If the true covering fraction of absorbing matter is as large as would be predicted if the dim 
continuum direction is correlated with the direction in which absorption is found, our picture of 
quasars changes substantially. The outward flow of matter moving at ~ 0.1c is no longer a small 
effect present in only a limited solid angle; it is, rather, a major element of the system, carrying 
substantially greater amounts of momentum and energy. Moreover, the velocity is no longer 
directed only within a small range of directions close to the surface of the accretion disk; there are 
fluid elements moving at very substantial angles relative to the disk surface. These considerations 
are likely to change our approach in dynamical studies of the BAL material. 

Beyond these dynamical issues, there are also a number of phenomenological consequences to 
this changed point of view. 

2.1. A new AGN X-ray population 

For only one BALQSO, PHL 5200, has an X-ray spectrum of good enough quality to 
determine an absorbing column been obtained; in this instance it is ~ 1 x 10 23 cm" 2 (Mathur et 
al. 1995). However, it has been possible to show statistically that most BALQSO 's are relatively 
weak in the soft X-ray band relative to their optical/ultraviolet flux. Interpreted as the result of 
absorption, the depression corresponds to columns ~ 10 22 cm -2 (Green & Mathur 1996). 

If the fraction of BALQSO's is nearer 0.5 than 0.1, and absorbing columns ~ 10 22 - 10 23 cm -2 
thick are a generic feature of BALQSO's, the number of quasars turned up in X-ray surveys 
sensitive to rest-frame photon energies greater than a few keV should increase by roughly a factor 
of two relative to the number found by soft X-ray surveys. Anisotropy in the X-ray emission is 
likely to be rather weaker than in the ultraviolet, but if present could modify this prediction. 
These new quasars should have hard X-ray spectra at rest-frame energies 1-5 keV (as a result of 
the absorption) , and most of them should also show BAL features in their rest-frame ultraviolet 
spectra. 
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Interestingly, the ASCA source count distribution shows exactly this effect: The number of 
sources seen by ASCA with fluxes ~ 1CT 13 - 1CT 12 erg cm -2 s" 1 in the 2-10 keV band is several 
times as great as the number that would be predicted on the basis of the corresponding Rosat 
counts if the sources had conventional quasar spectra (Cagnoni et al. 1997). Some of these may be 
AGN obscured by dusty tori as in Seyfert galaxy unification models (Madau et al. 1994; Comastri 
et al. 1995), but others may be BALQSO's. 



2.2. The radio-loudness of BALQSO's 

As already discussed by Goodrich (1997), a reduction in the optical flux of BALQSO's 
increases their ratio of radio to optical flux. He found that the (rough) statistics of known 
BALQSO's indicated that their optical flux had been depressed relative to the mean of all quasars 
by a factor of 3 - 5. 

More precisely, we wish to compute the ratio of radio flux (preferably extended radio flux, so 
that it is isotropic) to observed optical flux for quasars found in optically-selected samples, i.e. 
at fixed observed optical flux. If the radio flux F r = RF*, then for a single quasar, the ratio of 
observed radio flux to observed optical flux is F r /F t, s = R/ f(fx). Averaging over viewing angle at 
fixed F b s , the mean radio to optical ratio for BALQSOs is 

v _ jy- dn[R/f{pL)]d 2 N/dF obs d^ _ D jy- ^r~ 2 M 

{*r/K bs )BAL - ^ dfid 2 N/dFofjsdfI ~ dpfa-Hpy W 



so that 

(F r /F obs ) BAL _ 1 dflf a - 2 (p) 



(4) 



{F r /F obs ) all f BAL J 1 dflf a - 2 (fl) ' 

If either of the two models for /(/x) used in §1 is assumed with a value of that gives an observed 
BAL fraction of ~ 0.1, {F r /F obs )BAL is enhanced relative to the mean ratio for all quasars by ~ 2. 
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2.3. Scattering contributions to emission lines 

Many people have thought about possible contributions to the broad emission lines due to 
resonance scattering in BAL gas (Surdej & Hutsemekers 1987; Turnshek 1988; Hamann et al. 
1993; Hamann & Korista 1996). Indeed, some (e.g. Turnshek 1988 and Hamann et al. 1993) have 
tried to limit the covering fraction by constraints on these contributions. However, several special 
features arise in this model, which change the character of these arguments: the large covering 
fraction enhances the magnitude of the scattering contributions, although the limb-darkening 
partially counteracts this; and the width of the scattered line is diminished both by the preference 
for discovering face-on quasars, and by the relatively greater equivalent width of the intrinsic 
emission line as seen by the scattering gas. 

2.3.1. Emission line equivalent width 

If we consider an isolated resonance line (e.g. CIV 1549), the limb-darkening by itself can 
modulate the equivalent width of the intrinsic emission line by factors of several (as would similarly 
occur for isolated non-resonance lines such as H/3 or CIII] 1909). However, the combination of 
scattering and absorption diminishes the range of observed emission equivalent widths. Scattering 
increases the equivalent width of emission lines seen from the polar direction, which would 
otherwise be the direction in which the equivalent width was least, while absorption, of course, 
removes line flux from the equatorial direction. 

The quantitative character of these effects depends strongly on details of the geometry and 
kinematics. It matters whether the emission region is very small compared to the absorption 
region, or nearly fills it. The scattered line profile also depends on the velocity gradients within 
the scattering gas, and on the direction of the flow (helical outflow, as predicted by Murray et al. 
1995, has different properties from radial, for example). 

To qualitatively explore scattering effects, consider an extremely simplified model: Suppose 
that the emission line region is far inside the absorption region, the velocity of the absorbing gas 
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is purely radial, and the absorbing gas is black across some range of velocities j3 m in < P < Pmax, 
but transparent otherwise. With these simplifying assumptions, the scattered luminosity in the 
direction h is 

dL rPmax r+H* fi^ AT' 

-f=/ dp dp! \ d(f) ' g (n,h')-^5(P-P^, (5) 

dil J/3 min J-iL, Jo d\V 

where primed quantities refer to the intrinsic spectrum. Photons of initial frequency v 1 and 
direction h' are scattered with probability g(h, h') into direction h and frequency v by material 
with speed /?*, defined by the pair of equations 



v' = "oy/gfc (6) 
v , =v l=gj£* t (7) 

for rest-frame resonance frequency v Q . 

When /i* ~ 0.5, /3 m j n ~ 0, and (i ma x — 0.05 - 0.1, the observed emission equivalent width 
of an isolated resonance line seen in a non-BAL direction is typically enhanced by a few tens of 
percent relative to what it would be without scattering, and ranges over roughly a factor of two 
from the polar to the equatorial directions. Part of this enhancement in equivalent width is due to 
low broad wings. Depending on how the equivalent width is measured, in practise only a portion 
of these wings may actually be counted towards the emission line equivalent width, with the rest 
absorbed into the apparent continuum. Interestingly, if the maximum velocity in the absorbing gas 
is ~ 0.1c, the extent of the wings is comparable to the extent of the "very broad line" component 
identified by Brotherton et al. (1994). 

If attention is restricted to the red half of the emission line, this model predicts a greater 
contrast between BAL and non-BAL quasars. The emission equivalent width in BALQSOs 
should be some 50% stronger, on average, than similar lines in unabsorbed quasars. Weymann 
et al. (1991) found no strong differences between the emission-line properties of quasars with 
and without broad absorption; however, the dispersion of properties in each sample is large. 
An additional effect that might mask possible differences is that the intrinsic emission might be 
blueshifted relative to the systemic velocity of the quasar, so that non-BALQSO's appear to have 
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greater red-side flux. Investigation of such effects will have to await more sophisticated modelling 
of line profiles in limb darkened BALQSOs. 

2.3.2. Ratio of absorption to emission equivalent width 

The ratio of the absorption equivalent width to the emission equivalent width is sometimes 
taken as an indicator of the intrinsic covering fraction because scattering by itself would result in 
these being equal if the covering fraction were unity (e.g. Hamann et al. 1993). It is, of course, 
reduced by any true emission. The observed range of this ratio in BALQSO's is ~ 5 - 0.5, with a 
median value ~ 2 (Hamann et al. 1993). In the simple model defined in the preceding paragraph, 
this ratio also depends on both the intrinsic emission equivalent width and the range of velocities 
in the absorbing matter, decreasing for greater intrinsic equivalent width and increasing for 
increasing /3 max - For /i* ~ 0.5, it is 2 - 4 when /3 max — 0.1 and W\/X ~ 0.01 - 0.015. 

2.3.3. Emission line profiles 

Hamann et al. (1993) argued that if the covering fraction of scattering gas were large, the 
resulting emission line profile (in the non-BAL view) would be too broad. However, if the intrinsic 
emission line flux is originally isotropic while the continuum radiation is limb-darkened, this story 
changes. A wind in an equatorial wedge with a covering factor of 0.5 can intercept half of the 
isotropic line photons while intercepting a much smaller fraction of the limb-darkened continuum. 
The scattered C IV profile that results is then not much broader than the original line because 
the contribution from scattered line photons dominates that from scattered continuum photons. 
Scattering in the radial outflow can only reduce a photon's frequency in the rest frame of the 
photon source. 

We have explored this effect using the very simple model of scattering described above. 
Suppose, for example, that the absorption profile is generated by pure isotropic scattering, 
Pmax = 0.05, and ji* = 0.5. Let us further assume that the continuum limb-darkening is 
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given by ^(a 4 ), and that the intrinic C IV emission line, emitted by a point source at the 
center, is Gaussian with a dispersion of 0.01c and peak-to-continuum ratio of 1.5 when viewed 
along the pole. The resulting C IV profile, when viewed over a variety of non-BALQSO lines 
of sight, has a FWHM of 0.015c, a peak-to-continuum ratio ~ 2, and a shape parameter 
S = (FW 1 / imax — FWzi± may ) / FW H M ~ 1.0. All these numbers compare favorably with the 
C IV line-shape parameters measured by Wills et al. (1993) for quasars in the same redshift range 
as the bulk of the Weymann et al. (1991) sample. The red sides of the model C IV BALQSO 
profiles have similar shapes but somewhat higher peak-to-continuum ratios (~ 2.5) because the 
continuum is substantially limb darkened. 

2.34. Enhancement of NV 1240 

The pair of lines La and NV 1240 is a special case. Because the absorbing gas sees the 
intrinsic emission red-shifted, and the intrinsically strong La emission line is blueward of NV 1240 
by only ~ 6000 km s -1 , the flux available to scatter in the NV feature is substantially magnified 
relative to the pure continuum. For this reason, scattering can increase the NV 1240 emission 
equivalent width by factors of several (Surdej & Hutsemakers 1987), provided that the covering 
factor > 0.3 (Hamann et al. 1993). 

To quantify this effect, we again use the simple model already described. In this case, 
because NV 1240 and Lya are always blended, one cannot measure the equivalent width of the 
observed NV 1240 line by integrating across the profile after subtracting off a fixed continuum 
level. Instead, we defined the flux in each of the two lines by a least squares fit to a model of two 
Gaussians, one centered at 1240A, the other at 1216A. Specific numbers for the relative strength 
of the apparent NV emission due only to scattering depend on parameters, but only weakly. For 
= 0.5, limb-darkening law /2(m)> Pmin = 0, (3 m ax — 0.05 - 0.1, and La equivalent width from 40 
- 140A(a range taken from Baldwin et al. 1989), the scattered NV 1240 flux is 0.4 - 0.5 x the La 
flux. This ratio — with no contribution from intrinsic NV emission — is right in the middle of the 
range found by Baldwin et al. (1989) for their sample of quasars with z ~ 1.7 - 2.5. 
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Strong enhancement of the NV 1240 line is of particular interest in view of the fact that its 
flux, relative to both CIV 1549 and Hell 1640, increases on average by a factor of 2 - 3 comparing 
AGN at z ~ 0.1 with those at z ~ 2 - 4 (Hamann k Ferland 1992, 1993). If the incidence of 
BALQSO's also increases sharply with increasing redshift (as anecdotal evidence suggests), the 
scattering contribution may explain this systematic change in line ratios. Further support to this 
idea that the strength of NV 1240 is due to scattering in the BAL gas is given by the principal 
component analysis of quasar spectra by Francis et al. (1992), which showed that broad absorption 
is related to excess N V emission and broader C IV emission. 



3. Conclusions 

We have shown that if anisotropy in the continuum emission of quasars is correlated with the 
direction towards broad absorption line gas, optical flux- limited samples of quasars are severely 
biassed with respect to the discovery of broad absorption line quasars. If, as the simplest models 
suggest, quasar continua are limb-darkened at the selection wavelengths, BALQSO's are severely 
underrepresented in such samples. The inferred covering fraction of broad absorption gas could 
well be nearer ~ 0.5 than the usual estimate of ~ 0.1. 

Such large covering fractions of absorbing gas would have numerous consequences: Hard 
X-ray surveys may discover roughly twice as many quasars as would be predicted on the basis of 
soft X-ray surveys; the ASCA survey may already have done just this. The fact that BALQSO's 
are somewhat stronger radio sources than the average radio-quiet quasar would also be explained. 
In addition, resonance scattering in the broad absorption gas could contribute to the emission 
lines. NV 1240, because it is just 6000 km s _1 redward of Lya, would be strongly enhanced, 
providing an alternative explanation [i.e. other than elevated N abundance: Hamann & Ferland 
(1992, 1993)] for the strength of this line in high redshift quasars. 
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